Generalized Jaynes-Cummings model as a quantum search algorithm 
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We propose a continuous time quantum search algorithm using a generalization of the Jaynes- 
Cummings model. In this model the states of the atom are the elements among which the algorithm 
realizes the search, exciting resonances between the initial and the searched states. This algorithm 
behaves like Grover's algorithm; the optimal search time is proportional to the square root of the 
size of the search set and the probability to find the searched state oscillates periodically in time. 
In this frame, it is possible to reinterpret the usual Jaynes-Cummings model as a trivial case of the 
quantum search algorithm. 
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I. INTRODUCTION 

One of the most simple and interesting quantum mod- 
els that studies the interaction between radiation and 
matter is the Jaynes-Cummings model (JCM) [1]. The 
model considers the interaction between a single two-level 
atom with a single mode of the electromagnetic field. The 
coupling between the atom and the field is characterized 
by a Rabi frequency, and a loss of excitation in the atom 
appears as a gain in excitation of the field oscillator. The 
collapse and the eventual revival of the Rabi oscillation, 
described by the analytical solution of the JCM, shows a 
direct evidence of the quantum nature of radiation. The 
use of JCM has permitted to elucidate basic properties of 
quantum entanglement as well as some aspects of the re- 
lationship between classical and quantum physics. Since 
it was proposed, the pattern has been of permanent inter- 
est in the quantum theory of interactions. In the decade 
of the eighties it was discovered that the model exhibited 
highly non classic behavior, and the possibility of exper- 
imental realization appeared. The relative simplicity of 
the JCM and its extensions has drawn much attention 
in the physics community and recently in the field of the 
quantum computing [111] . In this work we use a gener- 
alization of the JCM to an N state atom interacting with 
a single field mode 

In 1994, Shor d described a quantum algorithm to 
decompose a number in its prime factors more efficiently 
than any classic algorithm. It was exponentially faster 
than the best known classical counterpart. In 2001 the 
experimental development of this algorithm has had a 
very interesting advance: Vandersypen et al. Q using a 
seven-qubit molecule manipulated with nuclear magnetic 
resonance techniques has reported the factorization of the 
number 15 into its prime factors 3 and 5. This algorithm 
illustrates a part of the theoretical challenge of quantum 



computation, i.e. to learn how to work with quantum 
properties to obtain more efficient algorithms. Tools like 
quantum parallelism, unitary transformations, amplifica- 
tion techniques, interference phenomena, quantum mea- 
surements, resonances, etc, must be used by the new com- 
putation science. Grover, in 1997, devised an algorithm 
[7J which can locate a marked item from an unsorted list 
of N items, in a number of steps proportional to y~N, 
that is quadratically faster than any classical algorithm 
Q. Continuous time search algorithms have been in- 
vestigated by a number of researchers [3, E(J EH ■ The 
essential content of these proposals is to built a Hamil- 
tonian (or alternatively a unitary operator) with the aim 
to change the initial average state \<j>) = ^2 m \(p m )/y/N 
into another state \ip s ), that also belongs to the same 
set {(v'm)} of N vector in the base of the Hilbert space 
0. This last state is recognized by an unitary operator 
called oracle, that is part of the global unitary evolution 
exp(—iHt), where the Hamiltonian is expressed as 



H=\4>){<f>\ + \ V .){<p a 



(1) 



The probability to obtain the searched state is 
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\{ip s \exp(—iHt)\(j)}\ , and it equals 1 after a time 7r-\A/V/2. 
In this frame, the search algorithm is seen as a rotation 
in the Bloch sphere from the initial average state to the 
searched state. Recently an alternative search algorithm 
was developed HI El El El that uses a Hamiltonian to 
produce a resonance between the initial and the searched 
states, having the same efficiency than the Grover algo- 
rithm. It can be implemented using any Hamiltonian 
with a discrete energy spectrum, and it was shown to 
be robust E3] when the energy of the searched state has 
some imprecision. The responses of this algorithm to 
an external monochromatic field and to the decoherences 
introduced through measurement processes was also an- 
alyzed in E3l- However we do not know of any experi- 
mental implementation, not even for a small search set. 
In this paper we present a resonant quantum search algo- 
rithm implemented with a generalization of the two-level 
JCM. 

The paper is organized as follows. In section |TT] we 
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consider a generalization of the JCM, in section IIIII we 
develop the search model. In Section |W] we present nu- 
merical results for our model. Finally in Section |V] we 
draw some conclusions. 



II. GENERALIZED JAYNES-CUMMINGS 
MODEL 

We shall consider the generalization of the JCM to an 
N state atom interacting with a single field mode with 
frequency lo synthesized by the following Hamiltonian [4] 



To calculate this transition probability independently of 
the initial and final numbers of photons the statistical 
weight of the photons must be incorporated, 

^ = ^EEp( n >H \{n Vl \W\m^\\ (10) 

n m 

where p(n) is the normalized photon number distribution 
and A is an unknown constant. Using Eq.© in Eq. (fTU|) 
we obtain the dependence of Pji with the average number 
of photons and the parameter Oo, 

p » = 4?{<™> + i + M™>}, (ii) 

N A 4 

H = hoj a t a+^£ /c ,S' fefe + -£!o ^ (a t Sjk + Skja) . (2) where (n) =J2 n P( n )- Taking into account the normal- 



ly 



k=l 



k=\ 



The photon creation and annihilation operators et^and a 
act on the photon number state \n) verifying 

[a,a f ] = 1, 

a'a\n) — n\n), 

a)\n) = \Jn + l\n + 1), 

a\n) — \fn\n — 1). 



(3) 
(4) 
(5) 
(6) 



£fc is the energy of atomic state \fk), is the atom-field 
coupling constant and it is fixed by physical considera- 
tions such as the cavity volume and the atomic dipole 
moment, Sjk is a transition operator acting on atomic 
states defined by 



ikmm) 



(J) 



where 5tm is the Kronecker delta. In what follows the 
subindex j shall indicate the initial state of the atom 
This state will be the starting state for the search 
algorithm and it can be chosen as the ground state for 
experimental purposes. 

Let us call \<p s ) the unknown searched state whose en- 
ergy e s is known. This knowledge is equivalent to "mark" 
the searched state in the Grover algorithm 0, [l(| EH- 
Our task is to find the eigenvector \(p s ) with transition 
energy io S j = {ej — e s ) /h from the given initial state \<fj). 
Then it is necessary to tune the frequency of the photon 
field with the frequency of the transition \(pj) — > \<Ps)- 
This means that the frequency of the cavity mode is se- 
lected as lo = u> s j. The transition between the atomic 
states is governed by the interaction term of the Hamil- 
tonian Eq.© 



W 



O N 

T ^ 

fe=i 



a) Sjk + S^a) 



(8) 



The transition probability to pass from the initial atomic 
state \tpj) with m photons to any final state \tpi) with n 
photons is proportional to Kn^lVFlmi/jj)! 2 . After some 
steps we get 

\{n(pi\W\mip.j)\ 2 = {(m + 1) 5 nm+1 + m £ nm _i . 

(9) 



N 



ization condition Pji = lj the dependence f2o with 

i=l 

the number of atom levels is obtained, 
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V(n) (N + l) + N v / (^' 



(12) 



The last step is valid for large N and (n). Note that 
flo depends on the photon distribution function only 
through the mean value of n. 



III. RESONANCES 

In the previous section we have determined the de- 
pendence of the atom-field coupling constant f2n with 
the number of atomic states N and the mean number 
of photons (n) . Now we want to study how this coupling 
constant determines the characteristic period of the dy- 
namics and subsequently the waiting time for the search 
algorithm. 

The dynamics of the system is given by the the 
Schrodinger equation for the wavefunction | 



ih 



d m)) 

dt 



H|*(t)), 



(13) 



where H is given by Eq.([2]). The global Hilbert space 
of the system is built as the tensor product of the spaces 
for the photons and the atom. Therefore atom-field wave- 
function is expressed as a linear combination of the 
basis {\tp m ) \n)} 



= ^2 X] bmn exp H ( £ m + wn) t/h] \<p m ) \n) , 

rn n 

(14) 

where {\n}} is the basis for the photons and {|v»ri)} is 
the eigenvector basis for the atomic Hamiltonian with- 
out electromagnetic field. The phase factor is introduced 
to simplify the final differential equations. Substituting 
Eq. (Ti"4"j) into Eq. (fT3")) and projecting on the state (<pi\ (k\ 
the following set of differential equations for the time de- 
pended amplitudes bik(t) are obtained 



Ik 



2i db, 
f2o dt 



= Vk exp [—i (u>ij 



t/h}b 



<jk-i 



3 



N 



5ij\/k + 1 22 b m k+i exp [-i (uj lm + u sj )t/H] . (15) 



This set of equations shall be solved numerically in the 
next section. Here we proceed to study their qualitative 
behavior. These equations have two time scales involved, 
a fast scale associated to the Bohr frequencies u>ik, and 
a slow scale associated to the amplitudes bi k (t). If we 
are interested in the slow scale all terms that have fast 
phase in the previous equations can be ignored; the most 
important terms are the ones with zero phase. In this 
approximation the previous set of equations Eq. (|15j) be- 
comes 



2? db sk 
Clo dt 



Vkbjk-i 



2i db 



f2o dt 



^ ~ VkTTb 



sk+1- 



(16) 



(17) 



These equations represent two oscillators that are cou- 
pled so that their population probabilities alternate in 
time. As we notice the coupling is established between 
the initial and the searched for states. To uncouple the 
previous equations we combine them to obtain 



df2 - + 4 b ]k , 



d 2 b 



O 2 

^-k^b sk . 
dt 2 4 



(18) 



(19) 



Solving these equations, for a given number k of photons 
and with the initial conditions bj k (0) — 1 and b sk (0) = 0, 
the following results for the amplitudes are obtained 



ci Q 



b jk (t) ~ cos( — Vfc + It), 



b s k+i(t) ^ -isin(-y"v/fe + It). 



(20) 



(21) 



Therefore the probabilities to obtain the initial and the 
searched states, independently of the initial number of 
photons, should be calculated as Pj(t) = ^p(k) \bj k (t)\ 2 , 

k 

Ps(t) — J2p(k) \b s k(t)\ 2 , which satisfy the conditions 

k 

Pj(0) = 1, P s {0) = 0. Averaging over the number of 
photons and using Eq. (fT2"]) these probabilities are 



(22) 



Pj(t) ~ cos 2 (Ct t) , 

P s (t) ~ sin 2 (CI t) , 
where the new angular frequency is 

CI- ^ 



N 



(23) 



(24) 



Then we see that the probabilities of the initial state 
Pj(t) and the searched state P s (t) oscillate harmonically 
with a frequency CI and period T = ^VN, while the 
probability of the other elements of the search set are 
neglected. If we let the system evolve during a time 



T 
1 



2A 



N. 



(25) 



and at this precise moment we make a measurement, we 
have probability 1 to obtain the searched state. It is 
important to indicate that this approach is valid only in 
the adiabatic approximation [l(| , this means that all the 
frequencies u> nm are much larger than Cl . Therefore the 
efficiency of our search algorithm is the same as that of 
the Grover algorithm and additionally it is independent 
of the number of photons. 

In the next section, we implement numerically Eq. (|15p 
and we show that this agrees with the above theoretical 
developments of the system in resonance. 



IV. NUMERICAL RESULTS 

The JCM has been used to understand the behavior of 
circular Rydeberg atoms, where the valence electron is 
confined near the classical Bohr orbit . This suggests 
to choose for our purpose an atomic model with an at- 
tractive potential whose quantum energy eigenvalues are 
e n = —eo/n 2 , where n is the principal quantum number 
and £o is a parameter. In this frame the Bohr transition 
frequency can be expressed as function of the parameter 
s a as u nm = s a (l/n 2 ~ l/m 2 ). 

We shall choose the numerical values of the parameters 
taking into account some previous experimental data. In 
Ref. [HI, [2(| the Rabi oscillation of circular Rydberg 
atoms was observed. The frequency of the single electro- 
magnetic mode was tuned with the transition between 
adjacent circular Rydberg states with principal quantum 
numbers 51 and 50, where the fast scale associated to 
this Bohr frequency was 1050,51 ~ 1007T GHz, that is very 
large compared with the fundamental Rabi frequency 
Clo ~ 507T KHz. In the micromaser configuration of Ref. 
[211 ] the field frequency of 21 GHz produced the Rydberg 
transitions used in this experiment, where the principal 
quantum number was about 63 and the Rabi frequency 
10 KHz. From the above we conclude that the ratio be- 
tween the Bohr transition energy and the vacuum Rabi 
frequency S = £o/^o to be used in our search algorithm 
should be taken as ~ 10 6 . 

We have integrated numerically Eq. (|15[) varying the 
parameter S in a range between 5 x 10 3 and 10 6 . The 
initial conditions are (i) a uniform distribution for the 
photons and (ii) bj k (0) is chosen in such a way that 
Pj(0) = 1. The calculations were performed using a 
standard fourth order Runge-Kutta algorithm. The 
procedure consisted in choosing at random the energy 
of the searched state and then to follow the dynamics of 
the set and of the initial state. We verified for several 
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FIG. 1: Probability distribution for all atomic levels at five 
different times. The search set {£}, with I — 1,2,...,N, has 
N = 50 elements (levels), 8 = 40000 and r is the optimal 
search time proportional to y~N. The initial state was taken 
to be j = f and the searched state s = 32. Note that the dis- 
tribution probability is essentially shared between the initial 
and the searched states. 



values of 6 that the most important coupling is between 
the initial and the searched states; the other couplings 
may be totally neglected. 

In FigQ]we show the probabilities for all levels at five 
different fractions of the time r. We see that the dis- 
persion among the states neighboring the searched state 
is relatively small for 5 — 40000. For higher values of 8 
the dispersion is even smaller, this confirms that our the- 
oretical approximation of two coupled modes is correct. 
Furthermore we conclude that the flux transfer process 
is essentially an interchange between the initial and the 
searched states and that the optimal time to measure the 
searched state is t. At other times we have less chance 
of perform a measurement of the searched state. 

Fig[2] shows the oscillation of the probability flux be- 
tween the initial and the searched states as a function 
of time for three values of S. The time is normalized 
for the theoretical characteristic time r. The evolution 
shows for the lower values of S an almost periodic behav- 
ior, however for the highest value 5 = 10 6 the behavior is 
completely harmonic. In this last case there is clearly a 
characteristic time when the probability of the searched 



FIG. 2: Probability distributions as functions of time for the 
initial (thin line) and the searched (thick line) states. The 
parameter S takes the values, from top to bottom, 5 x 10 3 , 
I0 4 and I0 6 . The size of the search set is N = 50. 

state is maximum and very near 1 and the initial state 
probability is near 0. The optimal time agrees with our 
theoretical prediction r. This periodic behavior and the 
proportionality between r and V7V are also found in the 
Grover algorithm 

V. CONCLUSIONS 

In this work we show how a generalized JCM to an 
N state atom interacting with a single field mode can 
be thought of as a quantum search algorithm that be- 
haves like the Grover algorithm; in particular the optimal 
search time is proportional to the square root of the size 
of the search set and the probability to find the searched 
state oscillates periodically in time. 

In the past, the biggest difficulty to build a JCM has 
been to obtain a single electromagnetic mode that inter- 
acts with the atomic transition. This difficulty has been 
overcome in the last decades thanks to the experimental 
advances in the handling of Rydberg atoms and to the 
building of micro-cavity for microwaves [U EH El, US 
l2lj . In the frame of this work we can interpret these 
devices as experimental realizations of the "analog" @ 
Grover algorithm in the trivial case of the search of a 
marked item in an unsorted list of 2 elements. However 
this new way of looking at the problem is different from 
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the usual point of view and opens new possibilities for the 
JCM. In summary, in this paper we reinterpret the JCM 
as the first step to build a more generic search algorithm 
with a large number of elements 
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